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ABSTRACT 



We analyse new spectra of the multiple system SZ Cam because previous studies found different values of the primary 
radial velocity amplitude. The older solutions of light curves also have different ratios of secondary to primary luminosity 
as inferred from the observed equivalent widths of spectral lines. We therefore reanalyse the light curves of the eclipsing 
pair. Only the light curve derived by Wesselink has a solution that agrees with the observed equivalent width ratio. 
The resulting parameters of the binary are discussed. Masses of Mi = 16.6 and M2 = 11.9 Mq, and radii Ri — 9.4 
and R2 = 5.4 R© are derived. We point out that radial velocities measured with the CCF method can be misleading 
when the method is applied to multiple systems with complex line blends. New radial velocities are also obtained for 
the visual component ADS 2984 A (HD 25639). 
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1. Introduction 



binaries: eclipsing - Stars: early-type - Stars: fundamental parameters - Stars: individual: SZ Cam, 



The early-type eclipsing binary SZ Cam (HD 25638) has 
been the subject of several studies in recent years. The bi- 
nary is a member of the open cluster NGC 1502 and is also 
' the northern component of the visual binary ADS 2984. 
With a combined spectral type of 9.5 and detached com- 
ponents, the binary is a potentially important source of ab- 
solute parameters of massive stars. However, the star was 
found to be a triple system - or, more specifica l ly, a s ystem 
consisting of two close binaries (|Maver et all 01)91;. The 
third contributor to the light makes the light curve solu- 
tion difficult, and also causes problems in the analysis of 
the complex spectra. 

The third body in SZ Cam manifests itself in various 
ways: as a third component in the spectral line profiles, by 
the presence of t he light time effect, a nd by additive light in 
the light curves (|Lorenz et al.lll998l hereafter LMD). The 
third body was also found in dependently by speckle inter- 
ferometry (iMason et al.lll998l). A new speckle measurement 
was published by IGorda et al.l (120071) . and the magnitude 
difference was given bv lBaleea et al.l lj2007f ) as m 95±0r 1 02 
at A535 nm. In the Washington Double Star CatalogQ, 
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SZ Cam is referred to as component Ea, and the speckle 
component as Eb. 

In addition to th ese studies, the s tar was analy s ed by 
[ Harries et all (11998L hereaf ter HHH). ICordal (1200(1 I2002L 
I2008D . and Micha lska et al.l (|2007h . HHH found some basic 
parameters of the system, namely the semiamplitude K\, 
which differs from those published by LMD. 

We collected new spectra with the 2 m telescope at the 
Ondfejov Observatory and also derived new solutions of 
the light curves. Based on these data, we expect to resolve 
the discrepancy between HHH and LMD. The already pub- 
lished results and our new results are listed in Table [TJ Note 
that the radial v elocity (RV) se miamp litudes K\ found by 
iMichalska et~aTl (|2007f > and by iGordal (|2008t) are close to 
those of LMD. In t his paper, the linear ephemeris given by 
IGorda et all (|2007ft is used: 

Primary minimum = HJD 2453676.1628 + 2 d 6984222. 



2. Spectroscopy 

The spectra analysed by LMD were obtained with the 
2.2 m telescope at the Calar Alto Astronomical Centre 
(hereafter CA). New spectra for SZ Cam were recorded 
with the coude spectrograph at the 2 m telescope of the 
Ondfejov Observatory. The region between 627 nm and 
678 nm is covered by 46 spectra with a two-pixel reso- 
lution of 12700 and a typical signal-to-noise ratio (S/N) 
of w 200. There are also two echelle spectra covering the 
region from 400 nm to 670 nm publicl y available at the 
ELODIE archive (|Moultaka et al.ll2004lh The spectral re- 
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Table 1. The published parameters of the RV curve (values in km s v ). 



Source 




m.e. 


K 2 


m.e. 


V 7 


m.e. 


Method 


Harries et al. 1998 


225.8 


3.8 


259.1 


4.0 


-12.0 


3.0 


CCF method 


Lorenz et al. 1998 


180.2 


2.0 


261.2 


3.8 


-2.9 


1.6 


MIDAS Gaussian fitting 


Michalska et al. 2007 


181.6 




268.2 








KOREL applied to Ha 


Gorda 2008 


192.0 


2.6 


266.4 


2.5 


-10.6 


2.0 


Gaussian fitting 


This paper 


189.4 


1.4 


264.1 


2.5 


-2.3 




Gaussian fitting 



1.00 



0.96 



0.92 



T3 
N 

E 

b 




1.00 



0.96 



0.92 



Phase 0.734 




300 



6670 



6680 
Wavelength (A) 



6690 



Fig. 1. Profiles of the line Hei 6678 A taken close to op- 
posite quadratures, at phases 0.285 and 0.734. The dashed 
lines represent the Gaussian fits of all three components of 
the system (the wide profile - primary, the shallow one - 
secondary, the narrow one - tertiary). The short vertical 
lines show the expected primary line positions according to 
the HHH solution. Four superimposed DIBs are also indi- 
cated. 



gion of the He I 6678 A line in two Ondrejov spectra taken 
at opposite quadrature phases is shown in Fig. [TJ The rel- 
atively narrow and strong third-body line is clearly visible 
in the centre of the blend. Owing to the strong reddening 
(Eb-v = 0.68), several DIBs are present: a stronger one at 
6672.0 A with equivalent width (EW) of about 0.02 A and 
several other weaker ones at 6681.1, 6684.4, and 6689.0 A 
(the wavelengths of these DIBs were measured in our spec- 
tra) Two short vertical lines mark the primary radial ve- 



2 At conjunctions, the DIBs 6672.0 A and 6689.0 A are not 
blended with the binary lines and can be measured relatively 
well. The other two DIBs are always superimposed on stellar 
profiles; their parameters were determined by optimizing the 
multiple Gaussian fit of the composite profile. 
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Fig. 2. Radial velocities of the eclipsing pair (the curves 
represent the solution with common Vy). 



locities (RVs) as they follow from the RV curve solution by 
HHH. These RVs are clearly incompatible with our spec- 
tra, and the large value of K\ published by HHH is not 
supported. 

We fitted the line profiles by Gaussians. The resulting 
RVs are listed in Table 2. Owing to the complicated struc- 
ture of the profile, reliable results can only be expected from 
spectra taken near quadratures, when the secondary line is 
separated. At phases around 0.75, the secondary compo- 
nent is blended with the DIB at 6672 A. Close to the Hei 
line, there is the Hell line 6683 A. We did not detect any 
contribution of this line, although in other stars of nearly 
identical spectral type - but probably of higher luminosity 
- the Hen line plays a role (e.g. at 6 Ori). 

We also refitted the CA spectra using an identical pro- 
cedure. The RVs of the eclipsing components are shown in 
Fig. The solution of the RV curves is listed in Table [3J 
Only RVs obtained using the line Hei 6678 A from spec- 
tra taken in Ondrejov and Calar Alto were used; velocities 
from phases closer than 0.08 to both conjunctions were not 
considered. In the covered spectral range, only the Ha line 
is present in addition to the 6678 A line. However, owing 
to its large width and disturbance by atmospheric lines it 
is not suitable for Gaussian fitting. 

In Table 2, RVs of the third lines are also contained. 
Since our new solution of the third-body orbit does not dif- 
fer considerably from that in LMD, we do not discuss these 
RVs here. Several spectra of the component ADS 2984 A 
were also obtained, and the RVs are given in Table 4. As 
the lines are wide, the resulting RVs are of low accuracy 
and the period of this probable binary has not yet been 
able to be determined. 



P. Mayer et al.: The O-type eclipsing binary SZ Cam (RN) 



Table 2. New radial velocities. 



HJD 




Phase 


RV 


RV 


Phase 


RV 


-2400000 


eel. 


pri. 


sec. 


third 


third 


53706 


2888 


0.1625 


-166.2 


216 


2 


o 


0757 


-18.5 


53706 


3075 


0.1694 


-164.9 


222 


1 


() 


0824 


-26.7 


53706 


3777 


0.1954 


-175.9 


249 


7 


o 


1074 


-21.2 


53706 


4290 


0.2145 


-191.4 


233 


3 


o 


1259 


-17.7 


53713 


1984 


0.7231 


198.4 


—261 


o 


o 


5460 


34.0 


53713 


2673 


0.7487 


179.2 


—261 


8 


o 


5707 


30.2 


53743 


2557 


0.8619 


152.3 


— 181 


g 


o 


2924 


14.8 


53744 


2390 


0.2263 


-182.5 


242 


3 


o 


6439 


21.4 


53747 


3102 


0.3645 


-140.4 


192 


g 


o 


7419 


-1.9 


53817 


4299 


0.3500 


-174.1 


197 


2 


o 


8118 


-20.9 


53818 


4625 


0.7325 


179.2 


—285 


2 


o 


1810 


0.5 


53832 


4657 


0.9220 


108.0 


— 131 


7 


o 


1875 


-0.7 


53938 


4909 


0.2135 


-182.6 


242 


7 


o 


0946 


-19.7 


53938 


5082 


0.2200 


-182.0 


237 


g 


o 


1007 


-18.1 


53938 


5465 


0.2341 


-182.4 


250 


4 


o 


1144 


-15.4 


53938 


5810 


0.2469 


-183.7 


262 


5 


o 


1267 


-15.0 


54014 


3879 


0.3400 


-163.3 


200 


5 


o 


2299 


-9.1 


54014 


4173 


0.3510 


-170.4 


203 


2 


o 


2404 


-6.4 


54020 


5757 


0.6331 


152.5 


—207 


1 


o 


4421 


25.7 


54028 


3257 


0.5051 


-42.8 


68 


6 


o 


2130 


-29.3 


54028 


3496 


0.5140 


-44.1 


68 


5 


o 


2216 


-25.6 


54028 


4000 


0.5327 


-21.7 


78 


4 


o 


2396 


-27.6 


54085 


4032 


0.6573 


163.9 


—216 


5 


o 


6199 


23.7 


54085 


4771 


0.6847 


178.2 


—241 


(i 


o 


6462 


25.4 


54085 


5942 


0.7281 


185.7 


—267 


2 


o 


6881 


15.6 


54085 


6617 


0.7532 


184.4 


—286 


1 


o 


7123 


12.2 


54085 


7320 


0.7792 


192.0 


—276 


2 


o 


7374 


12.2 


54097 


2111 


0.0332 


-51.2 


51 


5 


o 


8416 


-16.7 


54097 


2940 


0.0638 


-94.7 


87 


4 


o 


8712 


-17.2 


54097 


3726 


0.0930 


-107.9 


118 


8 


o 


8993 


-15.0 


54097 


4929 


0.1376 


-155.8 


213 


1 


o 


9423 


-28.5 


54097 


6107 


0.1813 


-180.6 


221 


6 


o 


9844 


-31.8 


54097 


7207 


0.2220 


-184.4 


253 





o 


0237 


-26.8 


54166 


4347 


0.6865 


188.5 


—286 


5 


o 


5909 


30.5 


54166 


4617 


0.6965 


188.8 


—269 


g 


o 


6006 


29.9 


54166 


5121 


0.7152 


186.4 


-277 


1 





6187 


27.3 


54166 


5434 


0.7268 


189.1 


-278 








6298 


27.3 


54533 


4269 


0.6890 


185.6 


-240 


3 





8011 


2.6 


54533 


4716 


U. /UOD 


lor (i 


-263 


6 





8171 


— 6.6 


54581 


3352 


0.4431 


-85.4 


109 








9297 


-15.5 


54590 


3348 


0.7785 


188.2 


-245 


7 





1473 


-2.0 


54594 


4061 


0.2874 


-182.4 


240 


9 





6029 


20.4 


54620 


3787 


0.9124 


114.8 


-142 


4 





8889 


-21.0 


54798 


6220 a 


0.9668 


35.8 


2 


6 





6160 


29.9 


54932 


4342 


0.5562 


63.2 


-94 


4 





4575 


22.1 


55060 


5181 


0.0222 


-55.2 


57 


3 





2512 


4.4 



Table 4. Radial velocities of ADS 2984 A. 



HJD 

-2400000 



53666. 
53670. 
53680. 
53797. 
53818. 
53938. 
53947. 
54028. 
54166. 



3657 

3298 

6565 a 

4386 

4075 

5452 

4385 

4483 

4878 



rev 


11.11 ) 


xtv 


„-l 
s 


-2400000 


km s — 


23.2 


54166.5829 


-31.1 


23.2 


54170.4173 


-8.4 


10.0 


54523.3472 


20.1 


12.4 


54593.3397 


18.3 


47.3 


54798.6220 b 


-17.2 


10.5 


55060.5444 


1.7 


13.3 


55074.4146 


-51.3 


26.0 


55075.6340 


-49.5 


15.6 







spectrum from the ELODIE archive 

both SZ Cam and ADS 2984 A exposed in same frame 



might be somewhat earlier than the type of the integral 
spectrum. Its luminosity is also higher, since the radius is 
larger than the radius of a MS member as shown in Sect. 
3; i ts log q is 3.7, a va lue typical of the luminosity class 
IV ([Martins et al.|[2005() . Although in stars of earlier spec- 
tral type the Hen lines are usually observed, this is not 
so in this case because the strengths of both the 6406 A 
and 6683 A lines stro ngly depend on logy. According to 
lLanz fe; Hubenvl (|2003h . the lines are absent for \ogg larger 
than 3.25 at T eS = 30500 K. 



3. New solutions of light curves 

In the following, we consider several SZ Cam light curves: 



— a photographic curve by IWesselinkl (|l94lh , the quality 
of which is comparable to photoelectric curves, since it 
is based on more t han 12000 measuremen ts; 

— UBV curves by iKitamura fc Y amasakJ (|l972t K-Y 
curves); 

— two curves obtained with int ermediat e wide filters cen- 
tred on 4720 and 5720 A bvlChocholl (|l980l ): 

— UBVR curves bv lGordal (pOOfT obtained with an area 
scanner, and published in a graphical form only. 



both SZ Cam and ADS 2984 A exposed in same frame 



Owing to the close vic inity of AD S 2984 A and E, 
Hipparcos photometry (|ESAlll99~7h is unhelpful. 

These light curves differ for two possible reasons: 

— the light curve is variable, or 

— the measurements are not reliable. 



their 



Among the existing classifications of S Z Cam, p r obabl y 
the most reliable is the one according to iBuddinsJ (|1978|h 
O 9.5 V. This i s the combined sp e ctral type. Following 
the scheme of IConti fc Alschulen (|l97lf h we confirmed 
this by comparing the lines 4471/4541 A (spectral type) 
and 4089/4143 A (luminosity) in the ELODIE spectra. 
However, a more detailed spectral type classification of 
the secondary and tertiary components using the ELODIE 
spectra is not possible due to weakness of the components' 
spectral features. These components are probably main- 
sequence stars, therefore their types and temperatures fol- 
low from their magnitude differences from the primary. 
Both should be later than the primary, so the primary 



The variability of stars of the given parameters is not com- 
mon. There are no indications that the magnitude in max- 
imum changes as it would be in the case of temperature or 
dimension changes. This explanation is therefore improba- 
ble. 

Observing SZ Cam with a standard PMT photometer 
was not easy, because of the close vicinity of the visual com- 
ponent. The K-Y and Chochol curves differ, each other and 
from the photographic curve. Gorda shows that there were 
deviations from the average curve in some nights. However, 
more data would be needed to verify the changes; the vari- 
ability of the comparison star might be also the reason. The 
photographic curve was collected over more than 6 years, 
which had the advantage of smoothing the eventual short- 
time light curve variability. 
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Table 3. Solution of the radial velocity curve of the eclipsing pair. 





Individual Vy 


Common V 7 


Parameter 


Primary 


Secondary 


Primary 


Secondary 


K (km s- 1 ) 


189.3±1.2 


263.4±1.8 


189.4±1.4 


264.1±2.5 


V 7 (km s^ 1 ) 


1.8±1.0 


-10.4 ± 1.6 


-2.3 


-2.3 


Average O — C (km s" 1 ) 


5.1 


11.5 


7.0 


12.6 


A ■ sini (Rq) 


10.09±0.09 


14.04±0.14 


10.10±0.10 


14.08±0.1 



Table 5. Solution of various light curves. 



Para- HHHT Chochol LMD^ Chochol Wesselink 

meter curves c curves' 2 



L£ 


27A%£ 


24.9% 


31.0%±1.5 


24.7%±3.8 


30.9%±1.8 


Q 


0.871 


0.871 


0.704±0.015 


0.710±0.064 


0.718 


Tx(K) 


29 725 


30 500 


30 500 


30 500 


30 500 


T 2 (K) 


27183 


27 200 


25 370±230^ 


27 270±650 


25 300±740 


i (deg) 


76.9±0.4 


76.7 


77.0±0.6~ 


77.0±1.8 


79.7±1.3 


L2/L1 in V 


0.773 / 


0.536 


0.422±0.031 


0.492±0.054 


0.262±0.029 


n/A 


0.319±0.010 


0.344 


0.380±0.030 


0.353±0.016 


0.379±0.004 


r 2 /A 


0.297±0.012 


0.270 


0.280±0.030 


0.265±0.021 


0.217±0.017 



a Solution by HHH; they used the two-filter data bv lChocholl (|1980I 1. 

LMD WD solution of Kitamura & Yamasaki UBV curves. 
c New solution of Chochol's curves with WD code and the mass-ratio q used by HHH. 

New solution of Chochol's curves with our q. 
e Assuming L\ + L2 + L-$ = 1. 

f Calculated from absolute visual magnitudes in HHH Table 6. 
9 Recalculated for T x = 30 500 K. 



A set of solutions is presented in Table [5] In our solu- 
tions Ti = 30500 K was assumed. A value close to this was 
already use d by HHH, and acco rding to newer temperature 
scales (e.g. iMassev et al.ll2005r i it is a reasonable value for 
the spectral type of the primary component. 

Naturally, as the curves differ, their solutions differ 
too. The largest differences are among the solutions with 
q = 0.71 and the HHH solution, with q = 0.871. We also 
solved the Chochol curves with q = 0.871. The adopted 
codes (LIGHT by HHH and WD by ourselves) are similar, 
hence the probable reason for the different results is that 
the Chochol curves constrain the system parameters only 
weakly. We also found that the scatter of normal points 
as published by Chochol is about twice as large as that 
of the UBV data used in LMD, and that the coverage of 
both minima is superior in the UBV data. When solving 
the Wesselink curve, we used the normal points given by 
Wesselink. The data and their WD fit are shown in Fig. [3] 

According to Table 3 in LMD, EW2/EW1 « 0.26 for 
the lines He 1 6678 and 4922 A. W hen EWs for these line s 
are calculated using models e.g. bv lLanz fc Hubenvl(l200l . 
it appears that they are practically constant in the given 
range of temperatures and are only weakly dependent on 
log g. Therefore, L2/L1 in V should be equal to the ob- 
served ratio EW2/EW1. 

Among the light-curve solutions, only the Wesselink 
curve gives such a luminosity ratio. Although we would 
gladly choose the solution of a photoelectric curve (and 
are not fully satisfied with the quality of our solution of 
the Wesselink curve because of the not too ideal fit of the 
minima bottoms) the spectroscopic evidence is clearly more 



consistent with the Wesselink curve, and is certainly incom- 
patible with the HHH solution. The corresponding param- 
eters of the system are listed in Table |6] 

Since the binary is well detached, we assumed that the 
evolution of the binary components was not affected by bi- 
nary interaction effects and can be well compared wi th evo- 
lution ary models for single stars. Using models by IClaretl 
(|2004h . with our values of R\ and T\ one gets the age 5.9 My 
and the mass 19.5 M Q . The age is close to 5.8 My given by 
iMalvsheval (|l997t based on uvby and H/3 photometry) for 
the cluster NGC 1502 and the mass is explainable by the 
often observed discrepancy between dynamical and evolu- 
tionary masses. In the case of the secondary, Claret's model 
gives the age 9 My and the mass 11.2 M Q for i? 2 and T 2 . 
In view of the possible errors, this model is acceptable. 
A larger R2, which would follow from other solutions in 
Table [5l would infer a yet older age. 

The mag n itude difference of the third body given by 
iBalega et alj (|2007l ) corresponds to the third body contri- 
bution of 29.4 %, i.e. agrees well with most of the solu- 
tions. Summing up the absolute visual magnitudes of all 
three components, one derives a total M v = —5.02. With 
^max = 6-98 an d A v — 2.27 (both numbers according to 
ICrawforc3ll994[ L the distance modulus is m — M — 9.73, 
close to Crawford's value of 9.71 (Crawford obtained this 
value as the average moduli of 20 cluster members using 
uvby and H/3 photometry analysis). 

It is certainly possible to assume synchronized rotation 
in this close binary. Accepting the solution of the Wesselink 
curve, the primary and secondary V TO t sini velocities should 
be 172 km s" 1 and 98 km s _1 . Using the Hei lines 4471 A 
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0.7 




-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Orbital phase 



Fig. 3. Solution of the Wesselink light curve with WD code. 
In the top panel, the differences "measurements minus the- 
oretical curve" are plotted. 



Table 6. Absolute parameters of the eclipsing system. 



Parameter 


Primary 


m.e. 


Secondary 


m.e. 


Semiaxis (Re) 


10.4 


0.1 


14.4 


0.1 


Mass (M ) 


16.6 


0.4 


11.9 


0.3 


Radius (R©) 


9.4 


0.2 


5.4 


0.2 


T cff (K) 


30 500 


fixed 


25 300 


740 


log L/Lq 


4.837 


0.005 


4.030 


0.009 


M bol 


-7.34 


0.04 


-5.33 


0.05 


BC (mag) 


-2.93^ 




-2.51^ 




M v 


-4.4.r 


0.04 


-2.8~ 


0.05 



a iMartins et all ((20051) 
b lLanz fc Hubeiivl l(2007l ). 



and 5876 A in the ELODIE spectra, we compared their 
FWHM with tables o flAbt et all ()2002l ) and a formula from 
M unari fc Tomasellal (11999). and derived V IO % sin i velocities 
of 160 km s _1 for the primary, 92 km s _1 for the secondary, 
and 60 km s _1 for the tertiary components, each with an 
accuracy of about 6%. The agreement is not perfect, but 
better than for other light-curve solutions. 



4. Conclusions 

We have attempted to explain the discrepancy between the 
K\ values of SZ Cam obtained by HHH and ourselves by 
the inadequacy of the CCF method for the solution of mul- 
tiple systems. HHH also found unrealistic velocities for the 
third body: absolute values of their third-body radial veloc- 
ities are larger than 4 km s~\ i.e. they disagree with the 
values given by LMD, ICordal (I2002L l2008h . and this paper. 
The CCF method produced erroneous results in the case of 
another binary in the HHH paper, IU Aur: no third light 
was found, although I3 « 30 % was deduced by various ways 
in this case. 

In a broader context, problems with t he application o f 
the CCF method were also highlighted bv lRucinskl (|2002l) . 
He mentioned namely the unreliable results achieved when 
the method is applied to three-body systems, and pointed 
out the uncertainty in the relative luminosity determina- 
tion - both problems now being illustrated in the case of 
SZ Cam. In the cases of SZ Cam and IU Aur, the low qual- 
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ity (S/N of only 40 and 20, respectively) and low number 
of spectra used by HHH might play a role. 

In several papers, the CCF method has been used to 
study extragalactic binaries. It is rather likely that a third 
component is often present in these binaries because of spa- 
tially unresolved field stars, and the results of the CCF 
method especially for binaries with light curves of small 
amplitude might be doubtful. 
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